Self-organization of atoms in a cavity field: threshold, bistability and scaling laws 



in 
o 
o 

< 

0\ 



J. K. Asboth 1 - 2 , P. Domokos 1 , H. Putsch 2 , and A. Vukics 1 
1 Research Institute of Solid State Physics and Optics, 
Hungarian Academy of Sciences 
H-1525 Budapest P.O. Box 49, Hungary 
2 Institute of Theoretical Physics, 
University of Innsbruck, Technikerstrasse 25, 
A-6020 Innsbruck, Austria 
(Dated: February 1, 2008) 

We present a detailed study of the spatial self-organization of laser-driven atoms in an optical 
cavity, an effect predicted on the basis of numerical simulations [P. Domokos and H. Ritsch, Phys. 
Rev. Lett. 89, 253003 (2002)] and observed experimentally [A. T. Black et al. in Phys. Rev. Lett. 
91, 203001 (2003)]. Above a threshold in the driving laser intensity, from a uniform distribution the 
atoms evolve into one of two stable patterns that produce superradiant scattering into the cavity. 
We derive analytic formulas for the threshold and critical exponent of this phase transition from a 
mean-field approach. Numerical simulations of the microscopic dynamics reveal that, on laboratory 
timescale, a hysteresis masks the mean-field behaviour. Simple physical arguments explain this 
phenomenon and provide analytical expressions for the observable threshold. Above a certain density 
of the atoms a limited number of "defects" appear in the organized phase, and influence the statistical 
properties of the system. The scaling of the cavity cooling mechanism and the phase space density 
with the atom number is also studied. 
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I. INTRODUCTION 



The manipulation of cold atoms and molecules by laser 
light is a rapidly growing field and has become a suitable 
ground for studying fundamental phenomena of physics 
both experimentally and theoretically 0. In the last 
decade, the emphasis has partly been shifted towards 
many-body effects in the dynamics of weakly interacting 
atoms 0, El • The mechanical action of the electromag- 
netic radiation held on free atoms rarely manifests these 
effects. The refractive index of a cloud of atoms is simply 
composed of the product of the single atom polarizability 
and the optical density. Standard laser cooling methods 
were also conceived on the basis of single-atom processes. 
Only at densities as high as the ones achievable in a 
magneto-optical trap does the dipole-dipole interaction 
between atom pairs give rise to a Lorentz-Lorenz-type 
refractive index and present an appreciable nonlinearity 
in the optical density 0,0- The underlying process, the 
reabsorption of spontaneously scattered photons in the 
atomic cloud [a, [3 , heats the atomic motion and hence 
limits the attainable minimum temperature. In addi- 
tion, this effect also introduces spatial instability into 
the atomic cloud and thus hinders degeneracy in phase 
space by optical means |8j . 

The mechanical effect of light on atoms inside a high- 
finesse resonator is substantially modified with respect to 
free space, which is the source of a variety of interesting 
phenomena in optical cavity quantum electrodynamics 
[3, EH • The basic reason is that a cavity photon makes 
many round-trips between the mirrors and thus the back 
action of the atom on the field, enhanced by the cavity 
finesse, cannot be neglected. As opposed to the external 
forces exerted by laser fields, the light forces in a cavity 
cannot be separated from the dynamics of the resonator 



mode, which is strongly influenced by that of the atom. 

The coupled atom-field dynamics can yield an efficient 
damping of the atomic motion via the mirror loss dissi- 
pation channel 0, 0] . Such "cavity cooling" schemes 
have recently been demonstrated experimentally |13lll4j . 
The fact that cavity cooling allows for replacing the spon- 
taneous emission, which is the dissipation channel in all 
laser cooling schemes, by irreversible photon loss from 
the cavity has important consequences. First, the in- 
ternal structure of the atom is not important and the 
mechanism can be operated on a wide range of species. 
Second, the problem of the reabsorption of spontaneously 
scattered photons, the source of the instability of atomic 
clouds at high densities, can be suppressed. 

The dynamics of atoms in a resonator is inherently a 
many-body problem even at a small density of the ensem- 
ble [H, ■ As all atoms are coupled to the same cavity 
mode, the modification of the field by one atom is expe- 
rienced by a remote atom as well as by itself. The cavity 
cooling mechanism may become inefficient since the del- 
icate dynamical correlation between one atom and the 
field mode could be perturbed by the motion of another 
atom [l7l llq . Indeed, one of the cavity cooling schemes 
was found to slow down linearly with increasing number 
of atoms . 

In a recent Letter we have predicted a cooperative be- 
havior of the atoms driven by a laser in a direction per- 
pendicular to the axis of a standing wave cavity po| . At 
high pump laser intensities (above a threshold) the ho- 
mogeneous atomic cloud self-organizes into one of two 
regular checkerboard patterns that maximize scattering 
into the cavity. The constructive interference of fields 
radiated by the individual atoms produces an intensity 
which depends quadratically on the number of atoms (su- 
perradiance) . Corresponding to the two patterns, there 
are two possible phases of the output field with 180 de- 
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grees difference, which have been observed in an experi- 
ment by Black, Chan, and Vuletic pi) . 

The onset of self-organization is relatively fast, on the 
microsecond time scale. A basic property of the present 
system is that the held created by the atoms traps and 
simultaneously cools them so that the organized pattern 
remains stable on a long time scale (10's of ms). The cav- 
ity cooling mechanism now acts on many atoms without 
losing efficiency. There is no external finite-temperature 
heat bath to define the temperature which, instead, is 
set by the dynamical equilibrium of the dipole force fluc- 
tuations and the cavity cooling effect. This is a distinc- 
tive feature with respect to the recently demonstrated 
collective atomic recoil laser in a ring cavity (CARL) 
HE Hdl HH, where a magneto-optical trap is nec- 
essary to stabilize the organized phase and the otherwise 
transient gain p(|, and also to inject noise for obtaining 
the phase transition-like behavior p?l I28L . 

In the present paper we discuss in detail the self- 
organization process, from the viewpoint of phase tran- 
sitions. A mean-field approach leads to a well-defined 
threshold in the pumping strength. Comparison to nu- 
merical simulations reveals effects scaling unusually with 
the atom density, a characteristic feature of this cavity- 
coupled many-atom system. 

The paper is organized as follows. In Sect.|nl the equa- 
tions of a semiclassical model are recapitulated, where 
the atoms are represented as simple linearly polarizable 
particles. Thereby the theory applies to a much wider 
class of particles than alkali atoms. The main features of 
the self-organization process, such as time scales, super- 
radiance, collective cooling, are surveyed using a numeri- 
cal example in Sect. Hill Then, in Sect. II VI we introduce 
a one-dimensional mean-field model and determine the 
threshold and critical exponent. In Sect. we present 
the results of detailed numerical simulations, which show 
effects beyond mean-field. The atom number enters the 
physics of the system in a form other than the density. 
Above a certain atom number, stable defects appear in 
the self-organized pattern and modify the system prop- 
erties, which is accounted for in Sect. IVII In Sect. IVIII 
the cooperative atomic behavior is discussed in detail by 
demonstrating the superradiance, and the ensuing im- 
provement of localization by collectivity. We conclude in 
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II. SEMICLASSICAL MODEL 

We consider N atoms in an open optical resonator 
(Fig. 2]). The atoms are illuminated from the side by 
a "transverse" standing-wave laser pump with frequency 
lu. This geometry corresponds to various experimental 
setups realizing the controlled transport of atoms from 
the side into a cavity 0, |3(j where the standing-wave 
pump amounts to a "conveyor belt" pll l3^ | . There is an 
efficient scattering of photons into the cavity mode quasi- 
resonant with the pump, wc ~ w, due to the enhanced 
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FIG. 1: The setup: transversely driven atoms moving in a 
standing-wave cavity. The laser intensity is given by the 
maximum free-space Rabi frequency (pumping strength) r\. 
The atoms couple to the cavity mode with 1-photon Rabi fre- 
quency g. The loss channels are: spontaneous emission (rate 
27) and cavity loss (rate 2k). 

dipole coupling described by the single-photon Rabi fre- 
quency g — Luy 2 [2eohV]~ 1 / 2 d eg , for a mode volume V 
and atomic transition dipole moment d eg along the cav- 
ity mode polarization. Large detuning of the laser from 
the atomic transition \lo — oja\ 3> 7, where 27 is the full 
atomic linewidth at half maximum, ensures that the up- 
per level of the atoms can be adiabatically eliminated. 
This model then describes a very general class of lin- 
early polarizable particles — in the following, we continue 
to use "atoms" for convenience. For the sake of sim- 
plicity, we restrict the atomic motion to two dimensions, 
along the pump laser and the cavity axis, coordinates x 
and z, respectively, without losing any relevant physical 
effect. Motion in the third dimension could be taken into 
account in the same way as along x. 

The quantum master equation for the density matrix 
reads 

P = ~[H,p]+jCp. (1) 
Here the Hamiltonian is 




N 



j=i i=i 

where a, a' are the boson operators of the cavity mode, 
Yj = (xj , Zj ) and Pj = (p X j , p Zj ) are the position and 
momentum vectors of the jth atom. The Liouville op- 
erator describing the cavity photon losses with rate 2k, 
and the spontaneous emission reads 

Cp = 2k ^apa 1 — —cJap — —pa^a^j 

N 

- To X)(^(ri)S(rj)p + p£ + (r,-)£(ri) 
3=1 

-2 J d 2 uJV(u)^(rj)e~ ttAUri pe ttxur *J5 t (r 3 -)) • (2b) 
In the above formulas, E(r) is the dimensionless electric 
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field, 

E(r) = f(r)a + rj(r)/g » cos(fcz)a + cos(kx)r] / g . (2c) 

The Rabi frequency of the driving laser is f]{r), whose 
position dependence is given by a cos(fcz) mode function 
for a standing-wave field. In the following, we are going 
to refer to the maximum value of the Rabi frequency r\ 
as "pumping strength" . The variation of the pump field 
along the cavity axis and that of the cavity mode function 
/(r) along the transverse direction (Gaussian envelope) 
are neglected. The detunings are defined as Ac = u — ujc 
and A a — oj — loa- The parameters 
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describe the dispersive and absorptive effects of the 
atoms, respectively, as they shift and broaden the res- 
onance line of the cavity. In the last term of Eq. \2h\. 
the integral represents the averaging over the angular dis- 
tribution N(u) of the random recoil due to spontaneous 
emission into the free-space modes. 

Instead of directly using the density matrix, we con- 
sider the evolution of the corresponding joint atom-field 
Wigner function |33j | . This can be systematically approx- 
imated by semiclassical equations for a set of classical 
stochastic variables, a, Pj, r^, the index j = 1,...,N 
labeling the atoms, 



relations, 

(£n£a) = ihT d n £ (rj)cos(fcZj) , 
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(UU) = 2h 2 k 2 T \£{r,)\ 2 ^5 nm + h 2 T 

d„S*(rj) d m £ (rj) + d n £(rj) d m £*(rj) 



(6a) 
(6b) 

(6c) 



where the indices n, m — Xj,Zj. The noise terms as- 
sociated with different atoms are not correlated. The 
complex dimcnsionless electric field £(r) is derived from 
Eq. H2c[). replacing the field operator a by the complex 
variable a. We iterate the coupled, stochastic Ito-type 
differential equations |0J by a Monte Carlo-type algo- 
rithm. 

There are two types of force terms in the equations 
of the momentum components. The terms in the first 
lines derive from the usual one-dimensional "optical lat- 
tice" potentials, the laser pump keeps the atoms inside 
the resonator via this term. In the second lines, the force 
terms originate from the coherent redistribution of pho- 
tons between the pump and the field mode. The po- 
tentials depending on the amplitude a, which itself is a 
variable, are not conservative (all but the optical lattice 
created by the transverse pump). The time-delayed cor- 
relations in the dynamics of the atomic motion and the 
field mode can result in a friction force on the atoms, 
known as cavity cooling [IJIMHIMIIllIllMlIi 



n =1 Ac— Uo cos 2 (fczj)J a— k+Tq cos 2 (fczj) 
j j 
— r/ c g cos(kzj) coa(kxj) + £ a , (4a 



Pxj = -hU Q (r]/g) 2 -^— cos 2 (kxj) 
d 

- ih(ii* s a - r] cfi a*)-g— cos(kxj) cos(kzj) + £ xj , (4b) 



Pz 3 = -hU Q \a\ 2 — cos 2 (fczj) 
ozj 

d 

- iH(r) eS a - ri cS a*)—cos(kXj) cos(kzj) + £ zj , (4c) 
ozj 

where the effective pumping strength for the cavity mode 
is 
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(5) 



These equations include Langevin noise terms £ Q , £ X j, 
and t^zj, defined by the non- vanishing second-order cor- 



III. SELF-ORGANIZATION 

We study he motion of the atoms in the cavity by nu- 
merically integrating the set of stochastic ordinary dif- 
ferential equations Q- To be specific, 85 Rb were consid- 
ered, with the 5 2 S!/ 2 ,F = 3 «-> 5 2 P 3 / 2 ,F = 4 transition. 
Starting from a gas of thermal atoms (random positions 
from a uniform, and velocities from a thermal distribu- 
tion) and no light in the cavity mode (a = 0), with the 
right choice of parameters we observe a buildup of the 
cavity field accompanied by the appearance of an orga- 
nized pattern in the spatial distribution of the atoms. 
This is illustrated in Fig. |2 where the trajectories of 
40 atoms during the initial 50 /is of a run are shown. 
The grid lines denote points of maximum coupling to the 
standing-wave cavity or pump field. Trapped atoms are 
oscillating about intersections of grid lines, where a sin- 
gle atom can scatter pump photons into the cavity mode 
most efficiently. For many atoms, however, destructive 
interference can inhibit the scattering process: the source 
term in Eq. I|4a|l contains the factor J^j cos(kxj) cos(fcz.,), 
which can be small even if all the atoms are maximally 
coupled due to the alternating signs of the summands. 
In contrast to this, in Fig. [21 only every second "max- 
imally coupled" site — the black or the white fields of a 
checkerboard — is occupied, leading to an efficient Bragg 
scattering of pump photons into the cavity. 



4 




-2-1 1 

cavity axis z [in units of X] 



FIG. 2: Two-dimensional trajectories of 40 rubidium atoms in 
a cavity, during the initial 50 usee. A checkerboard pattern of 
trapped atoms emerges, untrapped atoms move mainly along 
the cavity axis. Parameters: 7 = 20//isec, (g,n) = (2.5,0.5)7, 
atomic detuning Aa = — 5OO7, cavity detuning Ac = — K + 
NUo, and the pumping strength 77 = 5O7. 



The emergence of a checkerboard pattern of atoms 
with every second point of maximum coupling empty 
happens only due to the good choice of the parameters 
ensuring positive feedback, as explained in the following. 
Initially, in the random position distribution some atoms 
scatter into the cavity in a given phase and some with 
an opposite phase, and thus most of the scattered field 
is canceled. The dipole force (first term of Eq. JJEJ) at- 
tracts atoms towards antinodes of the pump (for red de- 
tuning, A a < 0), but almost no field in the cavity means 
no substantial modification of the uniform position distri- 
bution along the cavity axis. This can be seen in Fig. [2J 
where most atoms are well trapped along the transverse 
axis but some meander along the cavity axis. Due to sta- 
tistical fluctuations either the in-phase or opposite-phase 
scatterers will be in tiny majority, and a small cavity field 
does build up. The dipole force due to the cavity field, 
first term of Eq. l|4"c|) , now attracts atoms towards antin- 
odes of the cavity. The crucial point to consider is the 
interference of the cavity and pump fields, giving the sec- 
ond terms in Eqs. (0Jj,c). The product cos(kxj) cos(kzj) 
alternates sign between the black and white fields of the 
checkerboard. For a right choice of detuning Ac, there is 
a positive feedback and the atoms are attracted towards 
the "majority" sites and are repelled from the "minority" 
sites, due to the interference. 

The initial fast buildup continues over a longer 
timescale, with the kinetic energy of the oscillating and 
the free-flying untrapped atoms dissipated owing to the 
cavity cooling mechanism (for the transverse pumping 
case and for the chosen detuning Ac, it is the one de- 



scribed in Ref. 34]). This leads to an increase of the 
ratio of trapped atoms and to a stronger localization in 
the vicinity of the antinodes. Simultaneously, the co- 
herent scattering into the cavity improves, giving a slow 
increase in the cavity field intensity. The time evolution 
of the photon number is plotted in Fig. [5] for this self- 
organization process of 40 and 160 atoms. In the latter 
case, the photon number scale was rescaled by a factor of 
16. This way, the overlap of the two curves demonstrates 
the superradiance effect, i. e., the intensity is quadrati- 
cally proportional to the atom number. 




2 3 
time [ms] 

FIG. 3: The time evolution of the photon number in the cavity 
on a long time scale, for N = 40 and N = 160 atoms (note 
the different vertical scaling). The parameters are the same 
as in Fig. |5] 

In this phase transition- like process, the reduction of 
the kinetic energy is not a good characterization of the 
cooling efficiency. For the motion along the cavity axis, 
the transition from vanishing photon number to the "su- 
perradiated" light field is accompanied by a change of 
the heat capacity of the atomic ensemble since the num- 
ber of quadratic degrees of freedom changes from one 
(kinetic energy) to two (kinetic and potential energy). 
In deep harmonic traps, even a very low level of excita- 
tion can correspond to high kinetic energies. Thus even 
though the temperature can increase, cavity-induced dis- 
sipation increases the phase space density of the system 
by improving localization. An appropriate measure of 
this process is the effective phase space volume of the 
system (the inverse of the phase space density), mea- 
sured by the Heisenberg uncertainty product AxAp x /h 
and AzAp z /h for each degree of freedom. Here (Ax) 2 
and (Az) 2 are the averages of the squared distance of the 
atom from the nearest antinode along the direction x and 
z, respectively. For a harmonic potential, the dimension- 
less effective phase space volume amounts precisely to the 
mean number of excitation quanta. For untrapped atoms 
with mean kinetic energy (p 2 /m) = ksT, the phase space 
volume is Xy/rnksT / '(4\/37i) . 

In Fig. QJ the time evolution of the Heisenberg uncer- 
tainty product is shown for three cases. When the pump- 
ing strength is below threshold (N = 40, 77 = IO7), the 
uncertainty product is a constant. Here the spatial dis- 
tribution of the atomic ensemble remains uniform — the 
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transverse pump is too weak to induce any noticeable 
spatial modulation at temperatures kgT rs Hk. Thus 
(AzAp z ) = const, and (AxAp x ) = const, reveal that the 
temperature itself does not change in either direction. At 
the cavity cooling limit fc^T = Hk the numerical value 
for the phase space volume from the last paragraph is 
13.4ft, in accordance with the value along the transverse 
x direction. Along the cavity axis z, however, the phase 
space volume is above this estimate, indicating a higher 
temperature. Above threshold, the phase space volume 
transiently jumps to high values for both directions, and 
then it is gradually decreased. Compression is apparently 
more efficient along the cavity axis, here the phase space 
goes considerably below the value corresponding to the 
uniform distribution. In both plots the two curves corre- 
sponding to N = 40 and N = 160 are very similar, which 
manifests that the cooling rate of the ensemble is inde- 
pendent of the atom number. This is a very important 
observation, being at variance with the expectation that 
the efficiency of cavity cooling mechanism is reduced for 
increasing number of atoms. This prediction was made 
for a setup where the external pump field is injected di- 
rectly into the cavity. Apparently it does not apply to 
the transverse pumping case studied here. 

With the help of these few examples, we surveyed 
three important properties appearing in the dynamics of 
a laser-driven atomic ensemble coupled to a cavity mode: 
(i) the system rapidly self-organizes into a checkerboard 
pattern in a trapping field, (ii) which is generated by a 
collective, superradiant scattering into the cavity, and fi- 
nally, (iii) the energy of the atoms is dissipated at a rate 
independent of the number of atoms. This behaviour re- 
quires a sufficiently strong pumping strength, indicating 
the possibility of a well-defined threshold separating two 
different stability regions. This threshold is discussed in 
the next section within the framework of a mean field 
approximation. 
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FIG. 4: Evolution of the phase space volume, along the cavity 
axis (up) and along the transverse pump (bottom), on a long 
time scale for various settings, same parameters as in Fig. El 
The number of atoms is N=40 and N=160 (black line), and 
N=40 but pumped below threshold, 77 = IO7, for the dashed 
line. 



one-dimensional motion along the cavity. 



IV. MEAN-FIELD APPROXIMATION 

The essence of the self-organization process can be un- 
derstood on the ground of conservative mean-field forces 
acting on the atoms. This amounts to treating the cav- 
ity field as if it responded immediately to the positions of 
the atoms. Cavity cooling, which is directly related to the 
time lag of the cavity field, is absent in this model. More- 
over, the mean-field approach corresponds to the thermo- 
dynamic limit of the system: N — * 00, g — > 0, k = const, 
with Ng 2 = const. Physically, the limit can be thought of 
as taking larger cavities (cavity length / cav — > 00) filled 
with a gas of atoms of constant density (atom number 
N cx Z cav ). Due to the V~ x / 2 dependence of the cou- 
pling constant g on the mode volume V , one then has 
Ng 2 = const., neglecting the variation of the waist of the 
mode. Moreover, due to larger photon travel time be- 
tween the mirrors, the reflectivity has to scale like cx 
to keep k = const. For the sake of simplicity, we analyze 



A. Potentials 

Taking Xj = 0, j = 1, . . . , N, according to Equation 
(|4"cjl each atom moves in a potential 

V(z) = U 2 cos 2 (fcz) + U x cos(fcz) (7) 

composed of the sum of a A/2 periodic potential stem- 
ming from the cavity field and a A periodic one arising 
from the interference between cavity and pump fields. 
The potential depths are given by 

U 2 = N 2 (cos(fcz)) 2 hl U (8a) 
U x = 2N(cos(kz)) hl (A c - NU a (cos 2 (kz))) . (8b) 

These, in the mean-field approximation, depend on the 
position of the individual atoms only via the mean value 

1 N 

9 = (cos(kz)) = — ^cos(fczj) , (9) 
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which can be considered a spatial order parameter, and 
via the bunching parameter 



1 N 

B = (cos 2 (fcz)) = — y^cos 2 (fczj) 



(10) 



i=l 



The order parameter has characteristic values: (i) 
O rs describes the uniform distribution, (ii) — * ±1 
corresponds to a self-organized phase with atoms in the 
even or odd antinodes, respectively. Finally, Iq repre- 
sents the maximum number of photons each atom can 
scatter into the cavity: 



NToB} 2 + [A c - NU B]< 



(11) 



In equilibrium the spatial distribution of the atoms and 
the above averages are time-independent, which makes it 
possible to attribute a physical meaning to the potential 
(|7|). Since the potential depends on the position distri- 
bution, however, the system is highly nonlinear. 

For Uq < 0, obviously Ui < and the cavity field gives 
a potential with "even" wells at kz = 2nir and "odd" ones 
at kz = (2n + The interference term (oc cos(kz)) 
discriminates between the even and odd sites, raising the 
energy of one of them and lowering that of the other. If 
2 \U2 1 < \Ui\ this effect is so strong that V(z) yields a 
potential with wells at the even and hills at the odd sites 
- or the other way around, depending on the sign of U%. 

The sign of U\ is crucial. To simplify the dependence, 
let us require a cavity detuning Ac < —N\Uq\ so that 
the second factor in U\ is always negative regardless the 
momentary configuration of the atoms. To be specific, in 
the following we are going to use 



A, 



NU Q 



(12) 



If the atoms accumulate around the even (odd) antin- 
odes, then 9~+l ( w — 1 ) and the U± cos kz poten- 
tial is attractive at the even (odd) , while repulsive at the 
odd (even) sites. Therefore Eq. H12|) is the proper choice 
for positive feedback that makes the runaway solution of 
self-organization possible. 

Two more parameters describing spatial order are used 
later in this work: (i) the "defect ratio", the ratio of 
atoms closer to minority sites than majority sites; (ii) the 
"localization parameter" , the position variance (along 
the cavity axis and/or the transverse pump) 



V z = 



i N 

-Y 

N ^ \ 7T 



(13) 



where Zi is measured from the nearest antinode of the 
cavity mode function. A uniform distribution of atoms 
gives a defect ratio close to 50 % and a localization pa- 
rameter of 1/12. For a self-organized pattern, both pa- 
rameters approach 0. 



B. Canonical distribution 

We suppose that the phase-space distribution of the 
atoms factorizes to position and momentum depen- 
dence, the latter simply given by a thermal distribu- 
tion with mean energy fcsT. There is no external finite- 
temperature heat bath to set fcsT, it is instead deter- 
mined by the dynamics (@J through the equilibrium of the 
cavity cooling and the Langevin noise terms. This allows 
for a position- and time-dependent effective temperature, 
effects neglected in this model. For a far-detuned pump, 
cavity losses dominate spontaneous emission, and an es- 
timate kgT « Hk is provided by the Einstein relation. 
The spatial density of the atoms in the potential V(z) is 
then given by a canonical distribution, 

p(z) = ±cM~V(z)/(k B T)), (14) 

with the partition function Z = J exp(—V(z) / (kBT))dz 
ensuring that p(z) is normalized to unity. In our case the 
potential V(z) is a function of the density p(z), therefore 
this equation has to be solved in a self-consistent manner. 

A direct approach to solving (|14fl is to use it iteratively 
to determine the p(z) for given values of the physical pa- 
rameters. We set the temperature to the cavity cooling 
limit ksT — fin. Note that the temperature parame- 
ter T just rescales the pumping strength 77, e.g., taking 
ksT = 2Hk would correspond to increasing 77 by a factor 
of y/2. The results thus obtained after 100 iterations of 
(|14fl for an experimentally realistic example are shown in 
Fig. El (empty diamonds). There we plot the percentage 
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FIG. 5: Numerical iterations of the mean-field approximation. 
The percentage of atoms near odd sites after 10 (full circles) 
and 100 (empty diamonds) iterations is shown varying the 
pumping strength r\. The vertical line is at 77* of 11911 . The two 
insets show two steady-state position distribution functions at 
two different pumping strengths. Parameters: k = 7/2, A a = 
-5OO7, Ng 2 = 2OO7 2 Ac = —k — Ng 2 /\A A \, k B T = ha. 



of atoms near odd sites, i.e., with |(2m+l)A/2 — z\ < A/4 
for any integer ro, as a function of the pumping strength. 
Below a certain threshold in the pumping laser ampli- 
tude 77* (vertical dotted line) the atoms are distributed 
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uniformly, for stronger pumping this symmetry is bro- 
ken. Two examples of such self-organized position dis- 
tributions obtained by the iterations are shown in the 
insets. Note also that the convergence of the iterations 
is slow near the critical 77* (critical slowing down) , this is 
evidenced by plotting the results after 10 iterations (full 
circles) as well. 

The uniform distribution, p(z) = 1/A leads to = 0, 
B = 1/2, which give V(z) — 0, thus the distribution 
fulfills 114|) trivially for any values of the physical param- 
eters. To investigate its stability, we add an infinitesimal 
perturbation to it: 



,(")(*) = I + e5 (*)I, 



(15) 



with e infinitesimal, and the general A-periodic perturba- 
tion function 



rn — 1 



5> a 



(A m cos(mkz) + B m sin(mkz)) 



B m — 1 



(16) 



Since the spatial average of g(z) disappears, p(°'(z) is 
normalized to 1. Iterating (|14H once leads to the new 
density function 

pW (*) = \- Ne^^(Ac - yf/o) cos(fcz) + o(e 2 ). 

(17) 

To lowest order in e the only relevant perturbation is 
that proportional to cos kz. Stability requires that the 
first-order correction in e be self-contracting. Substitut- 
ing l|12f> for the cavity detuning, we have the following 
instability threshold for the uniform distribution: 



' N 

■V/,./' ( y \U Q \ 



K > kuT 



(18) 



For far-detuned atoms \Aa\ 3> Ng 2 /n, i.e., when the 
total cavity mode shift by the atoms is much smaller 
than the cavity linewidth N \Uq\ <C k, this translates to 
the following threshold on the pumping strength: 



77 > T] = 



k B T k \A A \ ^y 



Ng 



(19) 



This approximation of the critical value 77* is in good 
agreement with the simulations shown on Fig. giving 
77* = 35.47, which differs by less than 7 from the actual 
value. To make the physical content more transparent, 
this condition can be expressed in terms of the transverse 
pumping power density (energy/unit area/unit time) as 



Pin > k R T 



1 



4fc 3 



(20) 



3N/V ' 

As shown in the next Section, the condition 119|) amounts 
to requiring that the potential depth along the cavity 
axis of a self-organized checkerboard pattern exceed sig- 
nificantly the energy scale of thermal fluctuations. 



C. Critical exponent 

The stability analysis of the uniform distribution of 
atoms p(z) = 1/A has revealed the critical value of the 
pumping strength 77*. Moreover, we have seen that for 
77 « 77*, the relevant perturbation gives 



p(z) = 1/A + ecos(kz)/X. 



(21) 



Going beyond first-order perturbations, the above for- 
mula allows us to solve Eq. (f 1 4|) self-consistently. 

Pumping the atoms slightly above threshold, 77 = 
77* (1 + 8), with 5 <C 1, we expect <|2*T|) to give a good ap- 
proximation of p(z), with e depending nonlinearly on S. 
Substituting it into (|14|) we need to expand the exponen- 
tial to third order to obtain a solution to lowest order in 
the small parameters. This analytical calculation yields 
e oc (5 1 / 2 , i.e., above the critical value e increases from 
as the square root of the dimensionless excess pumping 
strength. The order parameter 9 = e/2 and the percent- 
age of majority atoms (2e/7r) both depend linearly on the 
small parameter e. Therefore, the analytical result shows 
the critical exponent of this phase transition to be h. 



V. NUMERICAL SIMULATIONS OF THE 
PHASE TRANSITION 

In the mean-field description of the steady state, the 
number of atoms enters only in the form of the atomic 
density N/V oc Ng 2 . The approximation is expected 
to be valid in the thermodynamic limit, i.e., N — > 00 
with the atomic density and the cavity loss rate con- 
stant. Trying to approach this limit in simulations of 
Eq. (0} we are in for a surprise. In Fig. El we show the 
measured percentage of defect atoms after 4 ms of simu- 
lation time as a function of the pumping laser strength. 
The thermodynamic limit is approached in three steps, 
N = 50, 200, and 800. The parameters are the same as 
in Fig. the atoms had random initial velocities from a 
thermal distribution with average kinetic energy tin. The 
initial positions were either uniformly distributed ( "up" ) 
or at "odd" points of maximal coupling ( "down" ) . In this 
controlled way we mimic the ramping of the laser power. 
Although the "down" curves show reasonable agreement 
with the mean field result, the "up" curves scale anoma- 
lously: a hysteresis is observed, whose breadth increases 
with the atom number. 

The hysteresis effect observed in Fig.[S]but absent from 
the mean-field prediction is due to the long timescales as- 
sociated with reaching a steady state. In fact, thermal 
fluctuations do not only alter the equilibrium solution by 
smoothing out the concentration differences due to the 
dynamics (this is correctly rendered by the mean field 
approach), but they can also delay significantly the on- 
set of that equilibrium. This slowing down is effective if 
the energy scale of thermal fluctuations is at least of the 
same order as the potential diffences due to the statistical 
fluctuations in the initial positions of the atoms. 
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Mean Field — 
N=50,g=2, up ♦ 
N=50,g=2, down -o- 
N=200,g=l, up » 
N=200,g=l, down -a- 
N=800,g=l/2, up ♦ 
N=800,g=l/2, down ❖ 



60 80 100 120 140 160 180 200 
T| [in units of y] 

FIG. 6: Ratio of defect atoms against pumping strength 77, 
4ms after the loading of the trap with a uniform ( "up" ) or or- 
ganized ( "down" ) gas of atoms. The different curves show the 
approach towards the thermodynamic limit. The parameters 
are the same as in Fig. 



The statistical fluctuations for a finite uniformly dis- 
tributed atomic ensemble lead to (N — 5N)/2 atoms 
around the kz = 2mr and (N + 5N)/2 atoms are 
around the kz = (2n + l)7r positions. Taking uniform 
distributions around the respective antinodes, we have 
5>os(fei) « 25N/tt and J2 cos2 ( kz i) ~ N / 2 - The Po- 
tential difference then reads 



AE = 2\Ui\ = hl 
. 4SN r] 2 g 2 



8SN 



N\U \/2) 



= h- 



7T KA\ 



(22) 



For the final expression in the second line, we consid- 
ered the far-detuned regime, \Aa\ 3> 7 ~ n ~ g, and 
N\U \,NT <£k. 

The self-organization occurs "instantly" if the trap 
depth originating from the statistical fluctuations ex- 
ceeds the average kinetic energy k B T of the atoms. Us- 
ing the expectation value of the finite-size fluctuations, 
SN w y/~N, we obtain 



V > m = 



k B T k \Aa \ \pH 
hn N 1 / 4 g 2 ' 



(23) 



Comparison with (|19(l gives 771 = y7r/8 iV 1 / 4 rf: the 
functional dependence of the two thresholds on the phys- 
ical parameters are the same, except for the "anomalous" 
scaling of 77-p with ./V as the thermodynamic limit is ap- 
proached. 

To check the laboratory timescale threshold for self- 
organization J22J), we performed simulations with the 
same physical parameters as in Figs. and starting 
from uniformly distributed atoms, but this time increas- 
ing the atom number with Ng 4 = const. In Fig. \7\ the 
percentage of defect atoms is plotted as a function of the 
pumping 77. These numerical results confirm that the 
threshold depends on Ng 4 , moreover, the value 771- = 837 
from l|23|) with k B T = Hk is also consistent with the sim- 
ulations. In the transition regime, there is a remarkable 



overlap of the curves corresponding to various number 
of atoms with keeping Ng 4 constant, indicating that the 
equilibration time also scales with Ng 4 . 



N=50, g=l. up 
N=200, g=0.7, up 
N=800, g=l/2, up 
N=3200,g=0.35,up 




60 80 100 120 140 160 180 200 
r\ [in units of y] 

FIG. 7: Ratio of defect atoms in the ensemble (measure of the 
spatial order) as a function of the control parameter 77 (given 
in units of 7). There is a clear threshold for self-organization 
which is independent of the atom number provided Ng 4 is 
kept constant. Parameters are the same as in Figs. HJ and 



The "down" curves of Fig.Elare invariant with respect 
to the thermodynamical limit, but the threshold they 
give is only 50-80 % of the mean-field prediction 77* of 
Eq. i|19|l • Again, this can be explained with an argument 
based on the comparison of energy scales. Instead of 
the statistical fluctuations in the position distribution, 
we now have a well defined initial state: ^ cos kzi = 
Y^cos 2 kzi = N. The threshold AE > k B T now gives 
the particularly simple result: 77 > r)i — t)*/2. This is 
slightly below the "down" curves of Fig. El the difference 
can be attributed to the nonoptimal coupling due to the 
position spread of the trapped atoms. 

In the last curve of Fig. N = 3200, the atom-cavity 
coupling g is the smallest parameter, and thus the sys- 
tem does not strictly realize the strong-coupling regime 
of cavity QED. In fact, self-organization does not de- 
pend on strong coupling: small g can be compensated 
by increasing the number of atoms. We remark that the 
experimental setup of was also operated out of the 
strong coupling regime. 



VI. STABLE DEFECT ATOMS 

The curve in Fig. [3 corresponding to N = 3200 devi- 
ates slightly from the other data and for large pumping 
strength it converges to a nonzero value of the defect 
atoms. In fact, for the parameters chosen there, de- 
fect atoms can be stably trapped at the minority sites 
of the checkerboard of maximally coupled points, which 
presents another important physical element of the sys- 
tem beyond the capabilities of the simple mean-field the- 
ory. In the following we discuss the condition for the 
appearance of defects and the upper limit on their num- 
ber. 
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For the analytical treatment of the defect atoms we 
use the ID model of Section II VI If the atoms are per- 
fectly self-organized, say, kzj — Irijix with integer rij 
for every j, then the potential depths of Eqs. I|8all scale 
with the number of atoms as U2 = —N 2 Io\Uo\, and 
Ui = -N2I K. Thus, for large enough N, the A/2- 
periodic potential is dominant and atoms can be trapped 
in the minority sites. These stable "defect" atoms reduce 
the superradiant scattering of the self-organized pattern. 

How many defects can stably reside in the pattern? For 
simplicity, we take N — M atoms exactly at kz = 2mr 
and M < N/2 "defect" atoms at kz = (2n + 1)tt, ne- 
glecting the position spread. We then obviously have 
J2 cos(kZi) = N - 2M and £ cos 2 (fcz,) = N. Substitut- 
ing this and the prescription Ijl2|l for the cavity detuning 
into (|8a|) we obtain 



U 2 
h 



-(N-2M) 2 \U Q \hI , 
-2(N - 2M)hI K , 



\VeS\ 



(24a) 
(24b) 

(24c) 



Defect atoms can persist if at every kz — nix there is a 
potential minimum: 

< 2 |[/ 2 | - |E7i| = 2hI {N - 2M) [{N - 2M) \U \ - k] , 



which entails 



N 
~2 



2 Un 



> M 



(25) 



If some defect atoms appear, the rise in their number is 
limited by the above inequality. In particular, if the left- 
hand-side is negative, there can be no stable defects: the 
condition for the possibility of stable defects reads 



N > N, 



thi 



\Uo\ 



(26) 



Note that this threshold is independent of the pumping 
strength. The maximum number of defects is limited by 



M < M„ 



N ~ N thl 



5 2 /|A, 



(27) 



we 



Working at large atomic detuning, \Uq 
find that defects can appear if the total mode frequency 
shift due to the atoms exceeds the cavity linewidth: 



Ng 2 >«|A A | 



(28) 



This inequality puts a lower bound on the atomic den- 
sity. Equivalently, it amounts to an upper bound on the 
atomic detuning A^: to avoid the occurrence of defects 
a large detuning can be chosen. This, |A^| 3> Ng 2 /n, is 
exactly the "far-detuned" limit of the previous Sections, 
used to derive the thresholds 77* , 77^ and r/i . Likewise, 
none of the curves in Fig. satisfy Eq. I|28|l . On Fig. [7| 



however, the curve corresponding to N = 3200 is above 
the critical density (|2SJ|. 

For comparison to the full solution of the dynamics, we 
numerically simulated the equations at fixed k = 7/2, 
g = 57/2, A^4 = — 5OO7, and 77 = 5O7. The number of 
atoms was varied from to 200, 25 runs with different 
random seeds were performed for each atom number for 
a duration of 5 ms. The conditions for self-organization 
derived in Sects. [^ an d[v] give for this parameter setting 
a threshold atom number N > 10-40 (for ksT = hn-hj). 
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FIG. 8: Defect atom ratio 5 ms after loading the trap as a 
function of the atom number. Each point is an average over 25 
runs of the simulation. The inset shows the average number 
defects. The physical parameters are: n — 7/2, g = 2.57, 
A A = -5OO7, f) = 5O7. 

The numerical results presented in Fig. |H1 show that 
the ratio of defect atoms averaged over the 25 trajecto- 
ries is well below 50%. For small atom numbers this is 
merely a "finite size" effect compatible with a balanced 
binomial distribution. For 15 atoms the defect ratio is 
still about half, for 18 atoms it is only a quarter of that 
expected from the uniform distribution. The ratio then 
drops down to the percent range, indicating stable self- 
organization, for 50 or more atoms. Concerning the ap- 
pearance of defect atoms, Eq. (|2l>|> gives N > N t hr ~ 40. 
Due to non-perfect bunching this threshold is shifted to 
somewhat higher values. The insert of Fig. |H1 shows the 
absolute number of defect atoms after 5 ms, averaged 
over the runs: the minimum at N ~ 60 followed by a rise 
can be identified with the threshold, which is in accor- 
dance with the previous, simple estimate. 

The transition from the perfectly ordered phase to the 
one where defect atoms can be present manifests itself in 
the statistical properties of the system. As discussed in 
Sect. IIIII the appropriate measure of the thermal excita- 
tions is the phase space volume, given by the Heisenberg 
uncertainty product AzAp z /h. In Fig. El this is plotted 
as a function of the atom number: both the final phase 
space volume of individual runs after 5 ms (dots) and 
the average over these trajectories after 2,3,4 and 5 ms 
(continuous lines) are shown. For very few atoms, the 
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FIG. 9: The normalized phase space volume of the system. 
Dots represent the values taken at individual trajectories at 
5 ms, the thick solid line is their average. The thin lines 
correspond to the average over the trajectories at 2, 3, and 4 
ms. The physical parameters are the same as in Fig. |H| 



phase space volume is scattered around the value of 13.4, 
corresponding to uniform spatial distribution and mean 
kinetic energy fcsT = Hk. For N as 50, the decrease 
of the average phase space volume and the reducing vari- 
ance around the mean show that the more stable the self- 
organized pattern, the less thermally excited the system 
is. At N as 60, a broad peak (due to the transition to 
a double well potential) heralds the appearance of stable 
defect atoms. This is followed by a slow but steady in- 
crease, which can be attributed to the rising percentage 
of defects. 

To indicate some of the dynamics, Fig. [§| shows the 
Heisenberg product at earlier times as well: at 2, 3, and 
4 milliseconds. The appearance of defect atoms slows 
down the equilibration process, but the respective curves 
converge uniformly to the one obtained at 5 ms. This 
demonstrates that the cooling time is closely independent 
of the atom number. 



VII. COLLECTIVE EFFECTS IN THE 
SELF-ORGANIZED PHASE 

The most direct evidence of cooperative action is su- 
perradiance into the resonator mode that can be mea- 
sured by detecting the power output from the cavity. In 
the self-organized checkerboard pattern each atom ra- 
diates with the same phase, and so the cavity photon 
number \a\ 2 increases quadratically with the number of 
atoms. This can be observed in the numerical simula- 
tions of Eq.(0}, gradually increasing the number of atoms 
as detailed in the previous Section. The cavity photon 
number is shown on a logarithmic scale in Fig. 1101 For 
many atoms in the cavity (N > 20) the intensity data 
are well approximated by the fitted quadratic function 
\a\ 2 = 0.08N 2 . The steady-state solution of Eq. iga} 
reads 
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FIG. 10: Photon number in the cavity 5ms after loading. At 
high atom numbers (TV > 20) the simulation results (dots) are 
fitted well by a quadratic function (black line), indicating that 
the atoms scatter cooperatively. The physical parameters are 
the same as in Fig. [H] 



If all atoms are perfectly at the antinodes 
((cos(kzj) cos(kxj)) 2 = = B = 1) this yields a 
coefficient of 0.125. The value from the simulations is 
30% below this: the difference can be attributed to the 
position spread and to the defect atoms. 

The superradiance has important effects on the spa- 
tial distribution of the atoms about their respective field 
antinodes. With increasing number of atoms the trap 
deepens so that the size of the atomic clouds about the 
antinodes is compressed, as shown in Fig. ^2 We note 
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(29) 



FIG. 11: The localization parameter T> z = J\ (kzi/%) 2 /N 
measured by the simulation (dots) and averaged over runs 
(continuous black line) shows an overall 1/N dependence as a 
function of the atom number. The dashed line is the approxi- 
mation of Eq. I3H . The physical parameters are the same as 
in Fig. |U 

again the appearance of the shoulder at N as 60 due to 
the transition into the double-well potential with stable 
defect atoms. Apart from this, the overall decrease of 
the localization parameter in the range N — 10, ... , 200 
is proportional to 1/N. As we show below, this scaling 
law can be derived by a careful examination of the self- 
generated trap potential. 
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We consider one-dimensional motion of the atoms in 
the limit where they are strongly localised in the vicin- 
ity of the antinodes. In this limit of harmonic oscillation 
the field amplitude in Eq. QJi) is coupled to the atomic 
positions only through the sum Y] zf, i.e., through the 
localization parameter defined in Eq. (|13|l . It is instruc- 
tive to introduce new coordinates in the configuration 
space: the mean radius r — zf /N, and a number of 
N — 1 angular coordinates tpj with canonically conjugate 
momenta p V j = Nmr 2 ipj, these latter decoupled from 
the field dynamics. Only the radial motion is damped by 
the cavity cooling mechanism, the angular ones are not. 
For the coordinate r the potential is composed of the har- 
monic attraction oc r 2 and a centrifugal repulsion oc 1/r 2 . 
There is a potential minimum and the radius is damped 
into it by cavity cooling. For large number of atoms the 
cloud size at an antinode can be simply estimated by the 
radius at the potential minimum. In this way we discard 
the role of fluctuations in this coordinate, assuming that 
the variance is much smaller than the mean. 

For a quasi-stationary field amplitude, the harmonic 
potential is mv 2 Nr 2 /2 with vibration frequency 

, = ^(!^M + W H'). (30) 

The angular kinetic energies are of the form 
p 2 0j /{2mNr 2 ), j = 1, . . . , (N - 1). The momenta 
p tpj can be estimated by their initial value in the 
unorganized phase, when all degrees of freedom have 
the same energy ksT/2 and the radius is r = \/i/48. 
The potential minimum is just at the radius where 
the harmonic potential energy equals the sum of the 
centrifugal energies. Simple algebra leads to the cloud 
size 

* ~£sJM& • (31) 

and since \a\ 2 oc N 2 , the 1/N law for the localization pa- 
rameter follows. Altogether, the increase in atom number 
results in a smaller cloud size in the vicinity of the antin- 
odes, a very important virtue of the collective atomic 
action. This compression is limited by the dipole-dipole 
interaction of the atoms near the same antinode, which 
effect was not taken into account in the present model. 

VIII. CONCLUSION 

A dilute cloud of non-interacting cold atoms in a high- 
Q cavity can undergo a phase transition if driven from 
the side by a laser sufficiently red detuned from an atomic 
resonance. Increasing the laser power above a threshold 
the atoms self-organize into a lattice so as to scatter most 
effectively into the cavity mode. In this way the atoms 
minimize their energy in the optical potential generated 
from the interference pattern of the cavity and pump 



field. Under proper conditions the atoms are cooled in 
this process giving long term stability to the pattern. 

The phenomenon has been previously seen in simu- 
lations and strong evidence was found in experiments 
plf. but threshold, scaling and efficiency of the effect 
with atom number, cavity parameters and system size 
remained unclear. Here we gave a thorough analytical 
discussion of the effect. A continuous density approach 
allowed us to derive analytical formulas for the critical 
pump power as a function of atom number and cav- 
ity volume and showed that the effect should persist if 
one scales up the volume at fixed atomic density, i.e., 
Ng 2 — const. Numerical simulations of the evolution for 
finite durations revealed a hysteresis between the ordered 
and homogeneous density phase on varying the control 
parameter, i.e., the pumping strength. This shifts the ob- 
servable threshold in the pumping strength, which then 
scales with Ng A . We still have a cooling mechanism for 
large numbers of any type of optically polarizable parti- 
cles. 

The system is composed of non-interacting atoms that 
communicate via a commonly coupled, single cavity 
mode. The cavity component of the system plays a mul- 
tiple role. First, the binding energy of the ordered phase 
is stored in the superradiantly enhanced field intensity of 
this single mode. Next, as an attractive feature of this 
system, the outcoupled field intensity directly serves as 
a possibility of time-resolved monitoring of the forma- 
tion of the ordered phase. Note that in setups without 
resonator, the uncontrolled scattered field can lead to a 
binding of micrometer-sized particles in an ordered pat- 
tern in liquid ^(JEI- In the cavity scheme, finally, the 
viscous motion is provided by the dynamically coupled 
single-mode cavity field. Apart from the geometry, this 
is also a distinctive feature with respect to CARL. 

From the laser cooling point of view, the picture is com- 
plicated by the phase transition. The appropriate mea- 
sure of the cooling efficiency is the phase space volume 
occupied by the system. The self-organization process 
reduces this volume below that of a system at the cavity 
cooling limit with uniform spatial distribution. More- 
over, the steady state is established within milliseconds, 
and this is closely independent of the number of atoms in 
the cloud (numerical simulations confirm this up to few 
hundreds of atoms). Collective cooling was previously 
known only for the stochastic cooling method 01 an d 
for common vibrational modes of trapped particles |43| . 
The collective behaviour strongly improves the localiza- 
tion, i.e., the size of the atom cloud pieces at the trapping 
sites is proportional to the inverse of the atom number. 

This work could be extended to various directions. 
First of all, the numerical simulations should confirm 
some of the statements of the present paper on a larger 
range of the atom number, or on a longer time scale 
(e.g., the dependence of the hysteresis on the duration 
of the evolution). Our prescription for the pump-cavity 
detuning Ac in Eq. I|12fl is probably impractical in the 
limit of large atom numbers, as one has to go closer to 
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the resonance in order to initiate the self-organization. 
Next, the maximum achiveable densities can not be de- 
termined in the present model as some of the limiting 
factors were omitted, e.g., vacum- mediated dipole-dipole 
interaction between the atoms, the eventual superradi- 
ant enhancement of the spontaneous scattering into lat- 
eral directions, etc. In extremely good cavities, quantum 
effects in the motion of the atoms begin to play an im- 
portant role which was not studied here. Finally, 
the nature of the phase transition is strongly determined 
by the geometry of the cavity mode: the possible trap- 



ping sites are defined by the antinodes of the cavity mode 
and the transverse pumping field. This situation can be 
essentially modified in a cavities with different geometry, 
e.g., ring or confocal resonator. 
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